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Abstract 

A simple sufficient condition on curved end of a straight cylinder is found that provides a local- 
ization of the principal eigenfunction of the mixed boundary value for the Laplace operator with the 
Dirichlet conditions on the lateral side. Namely, the eigenfunction concentrates in the vicinity of 
the ends and decays exponentially in the interior. Similar effects are observed in the Dirichlet and 
Neumann problems, too. 

Key words: thin domain, spectral problem, boundary layer, trapped modes, localization of eigen- 
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1 Introduction 

1.1 Formulation of the spectral problem 

Let uj C M n_1 be a domain bounded by a simple closed Lipschitz contour da; and n > 2. Let also H± be 
Lipschitz functions in Zu ~ lo U duo. Given a small parameter ft, > 0, we introduce the thin finite cylinder 
(Fig. P) 

fl h = {x = (y,z) S K' 1 - 1 xl : r] := h^y £ lu, ±z < 1 ± hH± (77)} (1.1) 
with the lateral side S h and the curved ends T^, 

S' 1 = dil h \ (WUrK) , T'l = {x:r 1 euj, ±z = l± hH± (77)} . (1.2) 
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Figure 1: Thin cylinder with distorted ends 



In the domain we consider the spectral mixed boundary value problem 

-A x u h (x) = \ h u h (x) , x&fl h , (1.3) 
u h (x) = 0, x e E h , (1.4) 
d n u h (x)=0, xeT'l, (1.5) 

where A x = A y + d 2 and A y are the Laplacians in R n and IR™ -1 , with d z = and <9„ stands for 
differentiation along the outward normal n defined almost everywhere on the Lipschitz surfaces Tjj_ . The 
variational formulation of the problem l]1.3 |) -(|1.5 p reads: to find A k eK and u h e H 1 (Q h ; S' 1 ) , u h ^ 0, 
such that the integral identity p] 

(V x u h ,V a: v) Qh =X(u h ,v) nh , v G H 1 (n h ;X h ) , (1.6) 

is valid. Here ( , ) nh is the natural inner product in the Lebesgue space L 2 (f2 ) and H 1 (yi h ; E 71 ) the 
Sobolev space of functions in fl h meeting the Dirichlet condition (|1.4p . 

The spectral boundary value problem p.3[) - p.5[) is of interest in many applicable disciplines related 
to acoustics and electromagnetism. Since h is a dimensionless parameter, the thin domain (jl.ip of length 
2 + O (h) can also be regarded as a long tubular domain with a transverse cross-section of unit diameter. 
In this way problem (|1.3|) - (|1.5|) can be understood as the spectral problem for a cylindrical waveguide 
with soft walls and hard ends (see [21 Ej and others). 

Since the left-hand side of Q1.6P serves as an inner product in H 1 (p, h ; Tr) and the embedding 
H 1 (p, h ] E' 1 ) C L 2 (£l h ) is compact, the spectral problem (|1.6p admits the positive unbounded sequence 
of eigenvalues 

< Ai < < \\ < ... < X'j < ... -> +oo (1.7) 

where the convention on repeated multiple eigenvalues is accepted. The corresponding cigenfunctions 
uf,U2) u 3i ■■■i u j > ■■■ m iP- h \ can be subject to the orthogonality and normalization conditions 

(u$,u$) ah =5 jlk , j, ke N = {1,2,3,...}, (1.8) 

where Sj k is the Kronecker symbol. According to the strong maximum principle, the first eigenvalue Aj" 
is simple while the corresponding cigenfunction may be chosen positive in CI . 

In the case H± = the cylinder has straight ends and the explicit dependence on h can be clarified: 

X^ q = h -^ p + ff, u^ q (x)=<p p (y)cos(^(z+lj). (1.9) 

The eigenpairs (|1.9p are renumerated with the two indices p, q G N. Furthermore, tpi, (p 2 , y>3, tfj, ... are 
cigenfunctions corresponding to the eigenvalues 

< hi < fi 2 < < ■■■ < fJ-j < ■■■ -> +oo (1-10) 

of the spectral Dirichlet problem on the cross section 

- \<P (v) = W (v) , V G 95 (y) = 0, y £ 9w, (1.11) 
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Figure 2: Thin curved trapezoid 



or, what is the same, 

(V^V^) w =M(^VO wl i/>&H 1 (u):=H 1 (u;du). (1.12) 
Similarly to (|1.8[) . the orthogonality and normalization conditions are satisfied: 

fe,^ = U j.fceN. (1.13) 

1.2 The localization of eigenfunctions 

All the eigenfunctions in (|1.9[) are oscillating and do not become infinitesimal as h — > +oo in any fragment 
LJh x (— i_,i_|_) of the cylinder = w/j x (— 1, 1) ; here u>h = {y '■ V € w } and 2± are arbitrarily fixed 
numbers, —1 < — Z_ < Z+ < 1. The main goal of the paper is to describe the localization effect for 
eigenfunction in the cylinder with distorted ends. In other words, we reveal the profile functions H± in 
(jl.ip such that at least the first eigenfunction vl[ concentrates in the vicinity of the ends and is of the 
exponential small order exp (— ft. _1 r) , r > 0, outside a neighborhood of Y±. 

A similar, but of other kind, localization effect appears in the Dirichlet problem for the Helmgoltz 
equation in a cylinder with the varying cross-section (see, e.g., [4]). To outline this kind of the localization, 
we briefly consider the Dirichlet problem in the thin curved trapezoid (Fig. [2]) 

n h = {x = (y, z) e M 2 : \z\ < 1, h^y e (0, H (z))} (1.14) 

where H E C 2 [—1, 1] is a positive function with the strict (b = —d 2 H (0) > 0) global maximum at z = 0, 
monotone for z € [— 1, 0] and z £ [0,1]. 

We present here only primary formal asymptotic analysis and refer to [SI [5] for details and much 
more scrupulous explanation and to (4j and [7] for different approaches based on the spectral theory 
of operators in Hilbert space and the r-convcrgcncc technique, respectively. Accepting the asymptotic 
ansatze for eigenpairs 

\ h = h- 2 ir 2 H{0y 2 + h- 1 A + \ h , (1.15) 
u h (x) = h~ 3 ^U (h~ 1/2 z) sin (h^nH (0) _1 y) + u h (x) , 

we use the standard procedure for the dimension reduction in the rapid variables r\ = h~ 1 y, ( = hr x l 2 z. 
We insert the ansatze (|1.15|) into the equation (|1.3|) and the Neumann boundary conditions (|1.4j) . Then 
we differentiate and apply the Taylor formula 

H { Z y 2 = H (Or 2 (l + z 2 H (0)^ b + O (|z| 3 )) = H (0)" 2 (l + hQ 2 H (0)" 1 b + O (tf |C| 3 )) 

(see, e.g., [FJ §4] and [5] for details). Factors on h~ 2 in the differential equation cancell each other due 
to the first formula (|1.15|) . Finally, factors on h^ 1 mould the limit ordinary differential equation 



- d 2 U (C) + B( 2 U (C) = AU (C) , C G M, 



(1.16) 



with the coefficient B = H (0) b > 0. Ei genpairs of the spectral problem (|1.16|1 . describing the harmonic 
oscillator, are known (see, e.g., [5]): 

Aj = B 1/2 (2j + 1) , (1.17) 
17, (C) = exp (^ 1/2 C 2 ) exp (-flV^) , j 6 R. 

The formulas (|1.17|l complete the asymptotic ansatze which exhibit the first asymptotic series of eigenval- 
ues in the thin curved strip (|1 . 14[) while estimates for the remainders \ h and u h in (| 1 . 1 5(1 can be derived 
by different approaches (see [U 03 [7] and others). Similar series of eigenvalues with stable asymptotics 
may be derived, e.g., by using local maxima of the function H in Q1.14|) . 

The distinguishing feature of the function Uj in (|1.17j) is the exponential decay o (exp (— <j( 2 )) for 
( — > ±oo; here a € (0, ^-B 1 / 2 ) . The corresponding eigenfunction u'j concentrates in the vicinity of the 
coordinate origin O and decays as o (exp (— /i _1 (tz 2 )) at a distance from O. On the other hand, the first 
eigenfunction of the problem (|1.3| - (|1.5jl in the rectangle (0, h) x (—1, 1) oscillates and is spread over 
the whole domain. In the trapezoid (|1.14|) with the concave upper side (Fig. 03) u\ is localized in the 
c/i 1 / 2 — neighborhood of the highest point in the upper side. 

Localization effects of this kind have been described in the papers 03 [4] 00 03, [10] for eigenfunctions 
of similar and other singular perturbed boundary value problems in domains in W 1 . 

In the present paper we study the localization effect of the other kind in the ch— neighborhood of the 
ends r j_ of the cylinder fl h . The localization of eigenfunctions relies upon the boundary layer phenomenon 
and the discrete spectrum of the limit boundary value problem (|2.3[) - (|2.5|) in the semi-infinite cylinders 

n± = = (t? ± ,C ± ) :»| ± 6w,( ± > -H± (?])} (1.18) 

obtained from the thin domain (|1.1[) by the coordinate dilation (|2.ip and the formal setting h = 0. It is 
known (see, e.g., [2]) can be readily verified that the continuous spectrum of the problem in Tl± coincides 
with the ray [fxi, +oo) but the discrete spectrum can appear below the cut-off [i\ which implies the first 
eigenvalue of the problem (|1.1|) . Points A* of the discrete spectrum in (0,/zi) give rise to the so-called 
trapped modes, i.e. solutions of the homogeneous problem (|2.3|) - (|2.5p with the exponential decay at 
infinity (see [11] and, e.g., reviews in 03 [T2] [13] ) . These trapped modes become the main asymptotic 
term in expansions of eigenfunctions of the problem (| 1 . 6|) . At the same time, the eigenvalues A* £ (0, fix) 
of the problem (|2.3p - (|2.5p in II± after the multiplication by h~ 2 approximate some eigenvalues in (jl.7p 
with the precision O (exp (— /i _1 r)) , r > 0. 

We especially mention the paper [5] (see also [5]) where both the approaches discussed above are 
combined to prove that the first eigenfunction of the mixed boundary value problem for the equation 
(|1.3p in a thin plate £l h C K 3 with the distorted lateral side may be localized near the point of the 
maximal curvature of the longitudinal cross-section. Note that in this case the decay rate is of order 
exp (— /i _1 r) inside the domain, but exp (— /i _1 / 2 cr) along the lateral side. 

1.3 Structure of the paper 

In Section [2] after a brief comment on the continuous and discrete spectra of the problem in the cylinder; 
we derive a simple sufficient conditions for the existence of eigenvalues in the interval (0,//i) below the 
continuous spectrum (see Theorems [2] and [3]) . 

In Section 03 we display asymptotic ansatze for eigenvalues and eigenfunctions of the problem (] 1 . 3jl - 
(|1.5p and estimate the asymptotic remainders (Theorem [4} . Our proof is based on certain estimates 
of Sobolev weighted norms of solutions to problems in Q h and II± (Lemmas 03 and 03) and differs from 
approaches in 03 El (Zl 01] ■ Also general results in the theory of self-adjoint operators in Hilbert space are 
used throughout the paper. 

In Section [4] we discuss similar spectral problems. First, we consider the two-dimensional problem 
(|1.3p - (|1.5p . for which the sufficient condition of nonempty discrete spectra only requires the negativity 
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Figure 3: The dumbbell domain 



of a coefficient in the Fourier series for the profile functions (0,1) 9 rj i— > H± (rj) (cf. formula (|4.2p ). 
We also find out a domain where the first eigenfunction concentrates at the both ends simultaneously. 
Then we consider the Dirichlet and Neumann problems in £l h for the equation (|1.3[) . A localization of 
cigenfunctions does not occur in the domain (| 1 . 1 1) with the Dirichlet conditions on the whole boundary 
(see Section |4~2")) . However, we demonstrate the same localization effect in the dumbbell domain (Fig. [3]). 
For the Neumann problem in the domain (jl.ip . only eigenfunction with large indices j = O (h J can 
admit the localization in the vicinity of r'_j_ and we detect the localization under a symmetry assumption 
on the domain Q . 

2 The boundary layer phenomenon 

2.1 The mixed boundary value problem in a semi-cylinder. 

We introduce the stretched coordinates 

H ± = (v,C ± ) = (h- 1 y,h- 1 (l T z)). (2.1) 

Since A x = h~ 2 A^±, the change of the spectral parameter 

X h i ► h~ 2 A (2.2) 

and the formal setting h — in the original problem (|1 .3[> - (|1.5[) leads formally to the limit problem in a 
semi-infinite cylinder (|1.18|1 with a curvilinear end 

-A£±U = AU in n±, (2.3) 

U = on S± = dIL±\TZ, (2.4) 

d,M = on T± = {e ± : r? G w, C ± = -#± (»?)} ■ (2.5) 

The variational formulation of the mixed boundary value problem reads: to find A g R + and U £ 
i? 1 (IT; S) , U 0, such that 

(Ve£/,VeV) n = A([/,TOn, ^^(n;S). (2.6) 

Here and in the sequel we omit the index ± in the notation. 

We relate the spectral problem (|2.3p ~ (|2.5p with the positive self-adjoint unbounded operator A in the 
Hilbcrt space L 2 (II) with the domain D (A) generated, according to [TH §10.1, 10.2], by the symmetric 
quadratic form 

a(U,V) = (V i U,V i V) n , U,V £ H 1 . (2.7) 

Remark 1 If the contour dw and the function H± are smooth, the operator A appears as the closure of 
the operator Aq with the differential expression -A( and the domain V (Aq) = H 2 (n)ni? 1 (II; S) . From 
a result in !15f based on the theory of elliptic problems in domains with piecewise smooth boundaries (see, 
e.g., f 1 6\ \ 1 7|/ and, particularly, U8\ \19 f ). shows that the domain T> (A) coincides with T> (Aq) in the case 
d n H± > on duo, i.e. for the acute edge j± = {£_ : rj S did , = —H± (?])} on the surface dH (cf. Fig. 
\Qa). If d n H± < and, therefore, the edge 7± is obtuse (cf. Fig. \Qb), then V (Aq) C V (A) . For any 
edge, T> (A) falls into the Kondratiev space U8f 



hi € H 1 (tt±;T ± ) : (1 + dist(£,j±)) k ~ 1 V|f7 G L 2 (II) , fc = 0,l,2} 
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Figure 4: Acute and obtuse edges 



2 N 2 N + 1 



2 N + 1 -1 2 N+1 



Figure 5: The plateau function 



2.2 The continuous spectrum 

Let us list properties of the spectrum a (A) of the operator A, which are well known (see, e.g., the review 
[3]) and readily follow from general results on semi-bounded self-adjoint operators in Hilbcrt space (see, 
e.g., [II])- The set a (A) belongs to the positive real semi-axis R+ in the complex plane C. The operator 
A has the continuous spectrum 

a c (A) = |>ti,+oo) (2.8) 

where [i\ > is the first eigenvalue of the Dirichlct problem (|1.12|) on the cross-section oj (see also 
formulas (| 1 . 1 1 [) and (jXTTUJ) ) . To verify this fact, one may observe that, for A > fix, the function 

(r,, C) » c-xp (±i (A - Ml ) 1/2 C) ipx (v) , (2.9) 

where tpi is the eigenfunction corresponding to /Xi and i the imaginary unit, satisfies the equation (|2.3[) 
in the entire cylinder ui x R and the Dirichlct conditions on the cylindrical surface doj x K but has no 
decay for both ( — > ±oo. Now multiplying (|2.9p by the plateau functions Xm with the graph in Fig. [5] 
provides the singular Weyl sequence for the operator A at the point A (see, e.g., [TH §9.1], [TH Thm. 
3.1.1], [3j [13] and others). This ensures that A belongs to the essential spectrum a e (A) . The kernel of 
the operator A — A is always finite dimensional (see, e.g., p~6]. Remark 3.1.5]) and, therefore, A £ a c (A) . 



2.3 The discrete spectrum 

The interval (0, Hi) C K+ contains the discrete spectrum ad (A) only. According to [21 Thm. 10.2.1], 
the lower bound of the spectrum a (A) can be computed as follows: 

min{A: A£a(A)}= inf a(U,U) 

u ei?Hn;£)\{0} \\U;L 2 {n)\\ z 

We emphasize that the infimum on the left is taken over all nontrivial functions in the domain H 1 (II; S) 
of the form (|2.7p which is bigger than the domain T> (A) of the operator A (cf. Remark [T]) . It clearly is 
one advantage of the theory p31 Ch. 10] applied here. 
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Our immediate objective becomes to find out a condition on the shape of the end T which provides 
a trial function W £ H 1 (II; E) such that 

a(W,W) </ii ||^;L 2 (II)|| 2 . (2.11) 

Indeed, the inequality (|2.11[) ensures that the lower bound (|2.10|) is less than the cut-off /ii for the 
continuous spectrum o~ c {A) and, hence, the discrete spectrum ad {A) cannot be empty. We employ an 
approach in [5j[20] (see also [13l §3.6]). For any e > 0, wc set 

W(£)=exp(-eOvifa)- (2-12) 
Since E = 911 n (du; x R) and ip x £ i? 1 (w; 9w) , the function W belongs to 1 (II; E) . Wc obtain 

/• p-\-oc 

||^;L 2 (n)|| 2 = / / exp(-2 £ C)^iW 2 dC^ = 

Jul J-H(rj) 

= 1 y exp (2eiJ fa)) ^ fa) 2 dry = (2.13) 
= ^+J H^cp, {rjfdri + e J H fa) 2 (pi fa) 2 dr, + O (e 2 ) , 

||V c W;L 2 (n)|| 2 = / / exp(-2 e C)(|V 1 ^ 1 ( ?? )| 2 + e 2 |^ 1 ( 7 y)| 2 )dC^ = 
= / cxp(2eH fa))(£i fa) 2 *? = 

J UJ 

= y £ + J H fa) |V^i fa)| 2 dr, + e(± + J^H fa) 2 |V^ fa) | 2 dr?) + O fa 2 ) . 

Here we have used the formulas 

||^i;L a (w)||=l, ||V^i;L 2 (w)|| = VMl, (2.14) 

which directly follow from ([Tig)) and fl~T2]) . 

Inserting (|2 . 13[) into (|2.1ip . we see that the terms of order e -1 cancel each other. Collecting terms 
of order e°, we observe that the inequality p. lip is valid with the function ()2.12j) and a small e > in 
the case 

J H( V )(\V ri < Pl (r ] )\ 2 -iJ 1 < Pl (v) 2 )dr 1 <0- (2.15) 

In other words, if one succeeds to find out a function H such that the relation (|2.15|) is met, one readily 
detects an eigenvalue Ai £ (0, /xi) of the operator A and simultaneously of the problem (|2.3p - (|2.5D . 

The function $ = |V^<pi| 2 — [L\ \<p\\ is of mean zero (see ([2ri4j)) and, therefore, changes sign in u>. 
Hence, the inequality (|2.15p surely can be satisfied by an appropriate choice of the profile function H. 
In this way H can be of rather arbitrary behavior (see Fig. [6j . 

Note that <& is positive in a neighborhood of the smooth contour duj because on one hand ipi = on 
dw and on the other hand |V,^i| > by the strong maximum principle. In Fig. 0a we locate a negative 
H in this neighborhood. The set {y £ u : <& fa) < 0} is not empty and includes the support of a positive 
H in Fig. [Sib. Finally, H is smooth and changes sign in Fig. [Sic. 

Theorem 2 If the function H £ C (Zu) meets the condition \2.15)) where {fi\,<pi} is the first eigenpair 
of the Dirichlet problem the mixed boundary value problem 12. 6\) in the semi-infinite cylinder 

n = {C = (ij I C):»?ew, C>-H( V )} (2.16) 

has an eigenvalue below the cut-off fj,i for the continuous spectrum (2.8\) . 
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Figure 6: Shapes which provide trapped modes 



2.4 A simplified sufficient condition for the existence of a trapped mode 

Let us simplify the integral in (|2.15|) under the smoothness assumption H £ C 2 (Zu) . Setting ip = Htpi 
in the integral identity (|1.12[) for the spectral pair {ixi, ipi} yields 

fii {(px, Hipx)^ = (V„<pi, V„ {Hpt))^ = (V^i.-ffV^i^ + {V v ipi, ^iV^-ff^ = (2.17) 

Thus, the condition (|2.15|) is equivalent to 

f ipi (rjf A V H fa) dr) < 0. (2.18) 

J UJ 

Theorem 3 // the function H £ C 2 (Zo) meets the condition \2.18}) where ipi is the first eigenfunction 
of the Dirichlet problem il.ll]) . the mixed boundary value problem \2. 3\) - ![2~~5\) in the cylinder i2.16\) has 
an eigenvalue below the cut-off [ii for the continuous spectrum i2.8\) . 

The inequality (|2 . 18[) becomes true for a subharmonic function H but false for supcrharmonic. Fig. 
[7J d and e, present two semi-cylinders (|1.18| determined with a subharmonic function H, while the 
condition (|2.15| is not so evident for application in these cases. At the same time, the semi-cylinder in 
Fig. O c, is given by a function H £ C 2 (to) which is not subharmonic but an eigenvalue A £ (0,/j.i) 
exists. 

If H is a harmonics, e.g., a linear function, the sufficient condition (|2.18[) is helpless. We attempt 
to make further use of the calculation (|2.13[) where we collect terms of order e. As a result, we see that 
(|2.11[) holds true with the function (|2.12p and a small e > provided 

\ + f H (n) 2 IV^x dr) — Hi f H ( v ) 2 ^ (r,) 2 dn < 0. (2.19) 

^ J UJ J UJ 

We now repeat the calculation (|2.17[) while changing H for H 2 and observe that the square of a harmonic 
function is supcrharmonic: 

l -A v H (t?) 2 = H (rj) A V H (t?) + \V V H ( V )\ 2 = \V V H ( v )\ 2 > 0. 

Hence, (|2 . 19|) converts into the false inequality 

|V„F(»7)|V ( V ) 2 d V <0 

and, therefore, a harmonic function H cannot assure the inequality (|2.11[) by the trial function (|2.12p . 



s 



Figure 7: The supcrharmonic and subharmonic ends of the semi-infinite cylinder 



3 Justification of the asymptotics 

3.1 The hypothesis and the theorem on asymptotics of eigenvalues 

We assume that the spectral problems (|2.3[) - ([2.5[) in the semi-infinite cylinders 1^ admit the eigenvalues 

A±...,A± ± g (o,/n-<y. (3.1) 

Theorems [2] and |3] provide sufficient conditions for the existence of the discrete spectrum below the 
cut-off fix and, accepting for H± either (|2.18| . or (|2.15|) . we fix Si > such that the inclusion (|3.ip is 
valid and N = N + U N~ > 1. In this section we prove the following theorem. 

Theorem 4 // the hypothesis (3. 1]) is true, the entries of the eigenvalue sequence jl.7]) of the problem 
(1.6]) in the thin cylinder (1.1]) satisfy the relations 



x'; - h- 2 A^ 



< c p exp (-Tph- 1 ) , h G (0, hp] , (3.2) 



where c p , t p and h p are certain positive numbers, p = 1, ...,N, and the eigenvalues (3.1]) are renumerated 
as follows: 

< A« < A (2) < ... < AW < tn -Si. (3.3) 

An information on the corresponding eigenfunction u\, ■■■,Uj^ is obtained in Section [3.51 as well. The 
eigenfunctions f/. of the problems (|2.6[) in Ii± can be subject to the normalization and orthogonality 
conditions 

(Uf,u£) n± =5 j , k . (3.4) 

After the renumeration of the eigenvalues A^ in (|3.3p . the corresponding eigenfunctions U^ p ) , of course, 
keep certain normalization and orthogonality conditions in n_|_ x II_ and we further still refer to (|3.4p . 
although a new way of writing the selfsame condition is accepted. 

3.2 A result on convergence 

We proceed with a primitive result on the convergence of the normalized eigenvalues in the lower range 
of the spectrum (|1.7j) and the corresponding eigenfunctions. 

Let u h G H 1 (fl h ; E' 1 ) be a solution of the problem ([0]) with 

A^G (0, h- 2 (m - r 2 )] , r >0. (3.5) 
Aiming to estimate Sobolev weighted norm of u , we introduce the positive weight function 

7? (<r )- / CXP ( T/l_1 (1 ~ |Z| ~ CliH ^ ' X G Cift\ 
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where t > 0, iltf = {x — (y, z) : 77 G u>, \z\ < 1 — cjjh} and cjj > max |if± (77)] is fixed such that fig C 
fi . The function (|3.6p is continuous and piecewise smooth while 

|V^r (a;) I < ThT l ll T (x) , a; G n'\ V X TZ T (x) = 0, ie!l\^. (3.7) 

We insert the test function w = 7Z 2 u G -ff 1 (fi' 1 ; S' 1 ) into the integral identity (|1.6[) and obtain 

A ft ||^ r ^;L 2 (O^) || 2 = A /l (u\li 2 T u h ) nh = (V x « h ,V x (ft£u fc )) n * = (3.8) 

= || (n T u h );L 2 (n h )\\ 2 - (u h v x n T ,v x (t^O)^ + 

+ (V x , l/Vx^r)^ - ||« fc Vx7Jr; L 2 f . 

Two terms on the right-hand side cancel each other. We integrate over z G (—1 + Cff/i, 1 — c#/i) the 
Friedrichs inequality in the small domain co = |y : G to observe that 

\\V y w;L 2 (n%)\\ 2 > h~ 2 /jn \\w;L 2 (fi&)|| 2 , w G H 1 (fl h ;E h ) . (3.9) 
Since K=loni£ f2' l \ft^ (see (pT6]l). from f[3?f |) - <[3T9 ]> it follows that 

A h ||w' 1 ; (O^) || 2 > A ft ||m' 1 ;L 2 (O fe \0^)|| 2 = X h \\ll T u h ;L 2 (Sl h \n%)\\ 2 = (3.10) 
= || V, (K T u h ) ; L 2 (il h ) f - ||u h V x ^ r ; L 2 (il h H ) f - A ft ||^ r u' 1 ; L 2 f > 
> h- 2 ( Ml - r 2 - /i 2 A' 1 ) ||ft T <A i 2 f . 

By virtue of the assumption (|3.5p , the factor on the last norm is bigger than h~ 2 (t 2 — r 2 ) and stays 
positive in the case r G (0,r ) . Thus, reading the relation (|3.10[) without its middle part, we see that 

\\K T u h ;L 2 (Ci h H )\\ 2 < h 2 X h (r 2 -r 2 )" 1 \\u h ;L 2 (il h )\\ 2 < (3.11) 
<to (^-r 2 )- 1 !!^;^ (O h )f. 

Returning back to (|3.8[) . we use (|3.1ip to conclude that 

|| Vx (K T u h ) ; L 2 (O fc ) f < (A' 1 + r 2 /^ 2 ) H^u*; L 2 (fi h ) f < (3.12) 
< h- 2 fa + r 2 ) (l + to (r 2 - r 2 )- 1 ) ||^; L 2 (fl*) || 2 . 

Taking into account that, by (|3.7[) and (|3.12D . 

||^ T V 2; ii' 1 ;i 2 (fi h )f < c(||Vx (K T u h );L 2 (n h )\\ 2 + h~ 2 \\K T u h ;L 2 (ft' l \^)|| 2 ) < ch~ 2 \\U T u h ;L 2 (fl h ) 
we formulate the obtained result. 

Lemma 5 Let {\ h ,u h } be an eigenpair of the problem U.6}) such that \\u h ; L 2 (Q h ) \\ = 1 and the 
inclusion VS. 5\) holds true. Then the following estimate is valid: 

h 2 \\R T V x u h -L 2 (tt h )\\ 2 + \\Tl T u h -L 2 (n h )\\ 2 < c(r), (3.13) 

where 1Z T is the weight function VS. 6}) with the parameter r G (0, ro) and the bound c(t) is independent 
of h G (0, 1] and {X h , u' 1 } but c (r) — ► +oo as t — > tq — 0. 
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Let the entry of the eigenvalue sequence (|1.7|) satisfy the condition ()3.5j) . From the corresponding 
cigenfunction u'j e H 1 (fl h ; we construct the two following functions in the semi-cylinders LI± 9 ^ : 

U-± (e ± ) = /^ 3/2 X (AC*) u* (to?, T /iC ± ± 1) • (3.14) 

Here ^ ± are the stretched coordinates (|2.ip and X G C°° is a cut-off function such that x(t) = 1 
for t < i and % (t) = for 4 > |, < x ^ 1- Note that ft 3 / 2 is a normalization factor caused by the 
coordinate dilation x >—>■ £ that together with Lemma [5] furnish the relation 

\\U^L 2 (U + )\\ 2 + \\U^;L 2 (n_)|| 2 = || x+ ^;i 2 (n")f + || x _^;L a (^)f = 

= i + ||(i-x+-x-)^;i 2 (^' l )H 2 > 

> l-ccxp(-(3/i) _1 r) ||ft r u);I/ > (f2 fe )f > 1 - Ccxp (- (3/i) _1 r 



Here x± (z) = x(lTz) and 1 — x + — x- on ly m ( — f > f ) ^ z where 1Z T (x) > c T exp \j3ti) rj , 
c T > 0, in accord with l|3.6[) . Thus, for a small ft, > 0, we have 

||L7jV;i 2 (n + )|| 2 + ||L/jL ;J L 2 (n_)|| 2 >i (3.15) 



By dSU), (j3~Tlj) and ifES]) . ([0]) . we obtain 



/i 2 Aj < Cj, 



I|c7-_f ± ; j? 1 (n±)|| ^ 3 y n (|v f (x^)r+|x<rj^= (3.16) 

= y ( (/^ \w x (x±4) 1 + ix±<r) < 

< c(ft 2 ||V. T ^ l ;L 2 (fi h )|| 2 + ||u);Z 2 (O ft )f) = c (ft 2 A^ + 1) < C r 
Thus, the following convergence occurs along an infinitesimal positive sequence {ftfe} feeN : 

ft 2 A^ 1 A° and U^ ± -> U° ± weakly in H 1 (H±; E±) . (3.17) 

Unfortunately, we cannot derive from (|3.17[) the strong convergence in L 2 (n±) in the unbounded domains 
However, we again make use of the weighted estimate (|3.13[) and write 

1 1 exp (tC^U^L 2 (n±)|| < Ijexp^ft-^lT^Xi^ 1 ^ 2 ^ h )\\ < 

<c\\TZ T u^;L 2 (n h )\\ < C. 

We now use the compact embedding H 1 (n± (T)) C L 2 (n± (T)) in the finite cylinder G n± : C ± < T} 
and the estimate 

\\U-±;L 2 (n±\n± (T))|| < cxp(-rT) ((exp^) u^l 2 (n ± \n ± (T))|| < 

< C* exp (-tT) 

with an infinitesimal bound as T — > +00. These yield 

U j± ~> strongly in i 2 (n±) . (3.18) 

We are in position to derive a variational problem for the tripple {A°, ?7j + , Uj-} ■ With any V± £ 
i? 1 (n±; E±) , we take the test function 

n h 3x^ v ± (x) = ft- 3 / 2 x± (z) V± (ft" 1 ?/, ft" 1 (1 T z)) 
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iii the integral identity (11.61) . Multiplying the identity with ft. 2 and going over to the stretched coordinates 
([2~1]) lead to 

= h 1 ' 2 (V a u£, V, (x±V±)) nh - h 1 ' 2 ^ ( x±u $, V±) a „ = (3.19) 

= (v*£tf ±l v^±) n± - /> 2 a£ (ttf ± , v ± ) n± + ft 1 / 2 (v^ ± v K x ± ) n „ - ft 1 /^ Kv,x ± , v,y ± ) n , . 

The first and second term on the right, by virtue of (|3.17p . (|3. 18|) . converge to (V|J7° ± , VgV±) n and 

A° (t/j±, Kt) n± i respectively. Since 7£(x) > exp ^(3ft) _1 T — c^^ for x Gsupp|V x x±l i the modulo of 
the third and fourth terms does not exceed 

cexp (- (3ft) _1 r) (ft 2 \\TZ T W ^ ; L 2 (fi ft ) f + \\Tl T uf, L 2 (fi h ) f) \\V±; H 1 (U ± )\\ . 

Hence, the limit passage k — > +oo, hk — > + converts (|3.19|) into the couple of integral identities (|2.6j) . 
(V 6 U9 ± ,V^V ± ) n± =A° j (u9 ± ,V ± ) n± , V± e H 1 (II±;E±) , 

while, in view of (|3 . 1 5|) . at least one of the functions ?7? ± is nontrivial. Since (|3.5[) provides the inequality 
A° < jiti — Tq , we conclude that A° coincides with one of eigenvalues in (|3.ip in the case Si > Tq . 
Let us formulate the obtained result. 

Proposition 6 // the eigenvalue X'j of the problem fl.6)) meets the condition i3.5\) . the convergence 
ft3.17\ ), (3.18\) occurs while at least one couple { A" , f/jj. } implies an eigenpair of problem S2.3\) - (2~5\) in 
the semi-infinite cylinder Tl± with the variational formulation i2.6}) . 

3.3 Approximation of eigenvalues and eigenfunctions 

Let H denote the Hilbcrt space H 1 (II±; E±) equipped with the scalar product 

(it, v) = (V x «, V x v) nh + hT 2 (u, v) ah . (3.20) 
We introduce the positive self-adjoint compact operator K, in TL by the formula 

(]Cu,v) = (u,v) uh (3-21) 

together with the new spectral parameter 

At = (A + ftT 2 )" 1 . (3.22) 

Comparing ()3.20|) - (|3.22p with ()1 .6|) . we see that the variational formulation of the problem (|1.3|l - l|1.5|l is 
equivalent to the abstract equation 

Ku = fiu, u £ H. (3.23) 

The operator K, has the essential spectrum {/j, = 0} (cf. [H] Thm. 9.2.1]) and the positive infinitesi- 
mal sequence = (A^ 1 + ft~ 2 ) 1 j of eigenvalues, while u 1 - = ||it^;H|| 1 Uj arc the corresponding 

normalized eigenfunctions. 

The following assertion is known as the lemma on "almost eigenvalues and eigenfunctions" (see [21] 
and, e.g., [14]). 

Lemma 7 Let (f> <E H. and a € K+ satisfy the conditions 

U;H\\ = 1, ||K#-tr#W|| = 9e (0,<r). (3.24) 
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Then K. has an eigenvalue in the segment [a — 9, a + 9] C R+. Moreover, for any 9\ € {9, a) , one finds 
coefficients a£, a l ^ +K _ 1 such that 

k+K-l 
j=k 

where ylj,, n\ +K _i is the complete list of eigenvalues of JC in the segment [a — 6\, a + 6i] C M+. 
The following couples ought to be chosen as approximate solutions of the spectral equation (13.211) : 

a h = ir 2 (Af + I)' 1 , (f> h = Wv^Hf 1 v^ ± , (3.26) 

where is an eigenvalue in (|3.ip with the corresponding eigenfunction UJ~ of problem (|2.6p in LI±, 
subject to the conditions (|3.4p , and 

v) ± (x) = h- 1 /\ ± (z) Uf (e-'y, e' 1 (1 T z)) ■ (3.27) 

To estimate the discrepancy 9 of the couple (|3.26p in the equation (|3.2ip . we need to study the decay 
properties of the eigenfunction Up . This can be made on base of the general theory of elliptic problems 
in domains with cylindrical outlets to infinity (see the papers [T51 [53] , ^ an< ^' e '^'' mon °g ra P ns 
|16l I17j). However, for the reader convenience we here present an elementary proof which is rather 
similar to our proof of Lemma [5] 

Lemma 8 Let {A = * = ,J7 =1= } be an eigenpair of the problem \2. 6\) in H± such that ||[/ ± ;L 2 (ILt) = 1 and 
the condition i3. 1)) is met. Then the inclusion exp (tC^) G H 1 (H±; S±) and the estimate 

||exp (rC*) E/* ; # ^Ili) || < c (m - A ± - t 2 )~ 1/2 (3.28) 

are valid with any r G ^0, — A^) 1 ^ 2 ^ . 

Proof. We omit the index ±. We introduce the weight function 

!exp(rc ff ), (<c H , 
exp(rC), (e[c H ,T], 
exp (tT) , C > T, 

which is continuous and bounded. Moreover, V x Rt,t = for £ ^ [cjj, T] and 

|V Sj R r , T (OI <TRr, T (t). (3.29) 

Inserting the test function V = R 2 T U £ H 1 (LI; S) into (|2.6j) and repeating the calculation ()3.8j) with 
an evident modification, we arrive at the relation 

A ||i? TiT [/;£ 2 (LI) || 2 = ||V c (f? T!T £/);L 2 (n)|| 2 - ||C/V ? i? r , T ; L 2 (n)|| 2 > (3.30) 
> Mi ||f?r,TE/;£ 2 (n ff )|| 2 - r 2 \\R T , T U;L 2 (U H )f , 

where II# = cu x (cjy, +oo) and c# = max{|if (77) | : 77 € U} . In (|3.30p we have applied the Friedrichs 
inequality in u> integrated over (ch, +00) 3 £ (cf. (|3.9|1 ) and the relation (|3.29|) . We finally obtain 

a\\u-,l 2 (n\n H )|| 2 > ( Ml -A-t 2 ) ||f? T , T c/;L 2 (n ff )|| 2 

or, by the normalization assumption on U, 

\\R T , T u- 7 L 2 (n H )\\ 2 < (^-A-T^Ap^L 2 (n\n ff )|| 2 < ( Ml - a-t 2 ) _1 a. 



n k+K-1 

E Kl 2 = 1 > ( 3 - 25 ) 

j=k 
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Since the function [c#,+oo) 3 ( m Rt,t (ViC) ' 1S monotone, the limit passage T — > +00 furnishes the 
inequality 

||cxp(rC)f/;L 2 (n ff )f < ( Ml - A - r 2 )" 1 A, 
which together with the relations 



\\R T ^(U;L 2 (n)f < c\\R TiT U;L 2 {n)\\ 2 , (3.31) 

Hexp^oc/;^ 1 (n\n H )f < \\u-h 1 (n\n ff )|| 2 < C (A + 1) , 

inherited from (|3~30l) and (JS^, lead to ([3~28]) . ■ 

To apply LemmaEl we need some calculations. First, using (13.28|) . we observe that 

\\vf;Hf = \J h X± (\V x Uf\ 2 + h- 2 \uf\ 2 ) dx+^J i W xX ± (2UfV x Uf + \uf\ 2 V xX ±) dx > 

(3.32) 

By the change of coordinates x 1— ► ^ ± , the first integral on the right reduces to 

5 X ± ( |v ^^' 2 + '^0 ^ L (1 - xl) (l v -^l 2 + ^ 2 l^l 2 ) dx - (3 - 33) 

The functions |V x x±| and 1 — x± vanish as ±z > | according to the definition of the cut-off functions 
X± and, therefore, bringing the weight cxp(r(' ± ) > cxp f(3/i) 1 in the last integrals in (|3.32p and 
(I3.33P yields the following upper bound for the integrals: 

cexp f (Shy 1 t)t f exp (rC ± ) f IV^C/fl 2 + ft~ 2 IC^H da; 

V ' ft J{xenh-.±z>2/3} v 1 J 1 / 

< cexp ((3ft)- 1 r) y exp (r^) (|VfC/± (£)| 2 + |C^= (£)f ) « < Cexp ((3ft)- 1 r 



The exponent r > is taken from Lemma [51 By ||J7- ; L 2 (n±)|| = 1 and (|2.6p . the first integral in (|3 . 33[) 



turns into | (l + Aj) . Thus, for a small ft, > 0, we have 



\vf;H\\>±(l + Af) 1/2 . (3.34) 



We further observe that 



= || £0 - a<j>; H\\ = sup I (/C0 - a<j); *) | = (3.35) 

-2 A ± f„.± 



11 ± ovir 1 

cr w ■ ; At sup 



where supremum is computed over all functions e such that ||^; 7i| = 1. Recalling (|3.27p and ()2.1D . 
we obtain 

V.*)^ - fc- a Af (^,*) n „ = h-^h- 1 {(V«E7± V C (x ± *)) n± - A± (^,X ± *) n± } + 

(3.36) 

+ ft- 1 / 2 (c/fv.xi,*)^ - ft~ 1/2 (v a t^, w„ x ±) nfc • 
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The expression in the curly brackets vanishes by virtue of the integral identity (|2.6|) with the test function 
£ i— ► x± i z ) * ( x ) which has a compact support and, therefore, falls into the function space H 1 (II±; £±) . 
Modulo of the last two terms in (|3.36[) does not exceed 

ch-^(\\U±;L 2 (wx [-§,§]) || ||V S *;L 2 (Q h )\\ + h\\V x U±;L 2 (w X [-§, f])|| hT 1 ||*;L 2 (fi' l )||) 

(3.37) 

< c/icxp (- (Shy 1 t) \\Uf\H 1 (n±)|| ||*;W|| . 



Thus, omitting the factor h which is unimportant since r < (/ii — A^) 1 ^"" is arbitrary we conclude the 
estimate 



< cexp 



(3.38) 



Owing to the relationship (|3.22| for the spectral parameters, Lemma [Jj extracts an eigenvalue A^ 
from the sequence (| X . T[) such that 

(\ h q +h-*y l -h 2 (Af + < cexp (- (3/1)" 1 r) . (3.39) 

This inequality transforms into the following one: 

| A* - h- 2 Af\ < ch-^cxp (- (3/1)" 1 r) (A* + l) (h 2 \ h q + l) < C7i- 4 cxp (- (Shy 1 r) (3.40) 

where the numbers C and t > depend on A^ only. Here we have used that A^ < fii (cf. (|3.1j) ) and 

A£ < h~ 2 Af + ch~ 4 cxp (- (Shy 1 r) (/ii + 1) (/i 2 A£ + l) =*> A£ < ch~ 2 (A^ + l) for a small h > 0. 

(3.41) 



We now assume that 



A(f-i) < A (p) 



(3.42) 



in the family (|3. 3|) . i.e., A' p ) is an eigenvalue with multiplicity x. Then the formulas (|3.26p with the 
eigenvalue A^ and the corresponding x eigenfunctions of the problems (|2.3p - (|2.5p in H± deliver linear 

independent approximate solutions m = p, ...,p + x — 1, for the spectral 

equation (|3.23|) . Furthermore, repeating the calculation Q3.32p and (|3.36p with obvious modifications 
yields the relations 



0« ) = S mik + O (cxp (- (Shy 1 r)) , m, ft = p, ...,p + x - 1. 



(3.43) 



We now set 0i = /30 in the second assertion of Lemma [7J where (3 > 1 is big but will be fixed 
independent of h. The condition 6i < h~ 2 (A^ + l) in Lemma [7J can be achieved by diminishing h. 
As a result, we obtain coefficients Oj such that 



fc+A'-l 

a l u l i ™ 

l=k 



k+K-l 



< 



where A£, Aj! +J< ._ 1 imply all eigenvalues of problem (|1 .€>[) such that 

(A ; ' l + /i- 2 ) _1 - /I" 2 (A'"' + <c/3cxp(-(3/i)" 1 r) 



(3.44) 



(3.45) 



We emphasize that, in comparison with Lemma [JJ we have enlarged in (|3.45p the bound 6\ = (39 
according to the previous estimate (|3.38p . However, this does not influence the conclusion: we chose null 
coefficients af m for new eigenvalues involved. 
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Recall that (vfr,v%) = 5 ltj . By (JSHJ) and J532J), wc have; 



k+K-1 

E ■ 

l=k 



Jim hr 




E a j' r; "i _ ^™> r ) - 



dm) Ar) 



(3.46) 



k+K-1 



k+K-1 



dm) 



E E « 




< cexp (- (3/i) V) +4/T 1 . 



Hence, for a small h > and a large /3 > 1, the coefficient columns a hp , ...,a hp+3< ~ 1 £ R k are "almost 
bi-orthogonal" that can happen in the case k > x only. As a result, we detect at least x eigenvalues in 
(|1.7|) subject to the estimate (|3.45[) . We fix an appropriate (3 and recall that, according to (|3.4ip . the 
estimate (|3.45p provides the inequality 



- /i" 2 A (p) ) < ch- A exp (- (Shy 1 t) < c p exp (-/i _1 t p ) 



(3.47) 



while Lemma [5] permits to take any r < (/ii — A^ p ') 



1/2 



with any r p < — and a certain c p . Note that we have replaced — by t p in order to subdue the factor h^ 4 

o 



Proposition 9 If the formula ^3.4%$ is valid with entries of the eigenvalue 113. 9\) . then at least x eigen- 
values \\ of the problem il.6\) satisfy the inequality \3.47\) - 



< 



3.4 The proof of Theorem [4] 

It remains to check up that, under the assumption (|3.42j) . the eigenvalues A p , A p+>t _ 1 and no other 
in (jTTTJ) meet the estimate p.47[) , i.e. k = p and K = x. We know that K > x. Assuming K > x, we 
infer the convergence h 2 \j — ► A^ for j = k, k + K — 1 while, by Proposition [H] and, in particular, 
formula (|3 . 18[) . the limits Uj± 6 H 1 (n±;£±) satisfy the problems (|2.6p in II ± and inherit the linear 



independence from u\ 



l k+K 



_ , . bincc A Cp) is an eigenvalue in 11+ of multiplicity x± and x = x + - 



the above inference and assumption arc invalid. In other words, K = x. 

If k > p, then, by Proposition [HI the total multiplicity of the spectrum (]1.7p in the segment 



c p exp 



(3.48) 



is bigger than p + x — 1 for a small h > and again we readily find out a contradiction with Proposition 
[6] Finally, the case k < p is impossible due to Proposition [9] which, dealing with all eigenvalues in (|3.3p . 
detect at least p + x — 1 eigenvalues Xf in the segment (|3.48|) . Thus, k = p and Theorem |4] is proved. 



3.5 Asymptotic expansions for eigenfunctions and the localization effect 

Let A(p) be an eigenvalue of multiplicity x (see (|3.42p ) in the list (|3.3p and let p + x — 1 < N. Then by 
Theorem 0] there exists h p > such that for /i £ (0, the segment 



(3.49) 



contains the eigenvalues Ap\ Ap+^i of the problem (|1.6p and is free of other entries in the eigenvalue 
sequence (|1.7p . Owing to the relationship (|3.22p , we see that the eigenvalues /z p = (A p + /i -2 ) 1 
(A p+ ^_ 1 + /i~ 2 ) 1 and no other eigenvalue of the operator /C fall into the segment 



2k 



+ A^^ + 2) \2/i 2 (A ( ^ +A^ 1) +2) 



(3.50) 
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We now apply the second assertion in Lemma [7] with B\ = h 2 6® where 0® is independent of h and, 
moreover, the segment [h 2 (A^ + l) - h 2 d\, h 2 (A^ + l) + h 2 6%] lies inside IJ53DJ). As a result we find 
coefficients af r such that, according to (|3.26p and (|3.34[) , 



p+*c—l 

E 



< 2— < c p /i" 2 exp (- (3/i) r) , l=p, ...,p + x- 1 



(3.51) 



We compare formulas (|3~20|l . (PT2"3"j) and |[LT5|) . ([0]) to derive that 



.ft. r2 



L 2 (n h )\\ =\\ l + h- 2 , i = p,..., P + h-i. 



Furthermore, by repeating the calculations (|3.32[) and (|3 . 36[) similarly to (|3.37[) . we obtain that 

Ivfrwf- fl + A«) < cexp (- (3/i)~M 



where v\ denotes the approximate solution (|3.27|1 corresponding to the eigenvalue A^) = A^ (sec (|3.26l) 



and cf. 



). Since A^ and Ajj, 



p+K—l 

E 

r—p 



1 satisfy (|3.39[) . the inequalities (|3.51[) convert into 



u n r -H 



< c p exp ^— (3h) 1 rj 



while, by means of the formula (|3.46[) with /3 = 0i/8, the estimates 



p+tc—l 

E 

q=p 



< c p h 



exp(-(3/i) r) , l,r=p,...,p + x-l, 



(3.52) 



(3.53) 



are valid. 

Now we employ the following simple algebraic fact (see, e.g., [25j Lemma 7.1.7], [26l Lemma 3.3]): 
under the condition (|3.53|) one can find out an orthogonal matrix b h = (b^ q ^J such that the inequalities 
(1 3 . 5 2 [) ensure the estimates 



- h 



p+K—l 

; E 

r—p 



< Cph 2 exp 



(- (3/i) 1 t \ , I =p,...,p + 



x-l. 



(3.54) 



Theorem 10 // A^ p ' is an eigenvalue of multiplicity x (see 13. 3\) and \3.42ty ), then there exists the 
coefficient columns b\ L = \ b^ v ', 6j ip+><_1 ^ , I = p, ...,p + x— 1, composing an orthogonal matrix of size 

are eigenfunctions of problem ll.6\) 



k x k and furnishing the inequalities \3. 54\ ) where 
under the normalization and orthogonality conditions 111. 8]) and v p , Vp +>c _i are the functions \3.21ty 
constructed from eigenfunctions of problem \2. 6\) in the semi-infinite cylinders Tl± (see (1.25)) under 
the normalization and orthogonality conditions {3.J$. 



Since the eigenfunctions Uj- decay exponentially in the semi-cylinders H± as £ 



[5]), the functions v' p l 



U 3± 



in (|3.54p and (|3.27p are of order h~ 



+oo (see Lemma 
1 / 2 in the vicinity of the end of the 



thin cylinder £l h but become exponentially small at a distance from this end. Thus, the estimate (|3.54[) 
obtained in Theorem [TU1 exhibits the localization effect discussed in Section [L~2l On the other hand, the 
structure of the weight function (|3.6p . which is exponentially large in the middle of the thin cylinder 
([l.ip , and the estimate (|3 . 13[) of weighted norms of an eigenfunction u h of problem (|1.6p ensure the 
same effect. We emphasize that in the case when A(p) = Af is a simple eigenvalue in the list (|3.3p . 
Theorem [10] provides the localization of the corresponding eigenfunction Up in a neighborhood of only 
one end while Lemma [5] cannot distinguish between the ends and (sec (|1.2p ). In Section |4~T1 
we demonstrate an example of the principal eigenfunction u\ which does not become exponentially small 
near the both ends. 
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4 Similar spectral problems 

4.1 Splitting of a multiple eigenvalue in the case n = 2 

If uj = (0, 1) C M 1 and 

il h = {x = (y, z) e M 2 : r) = h^y € (0, 1) , z < (-1 - hH_ (rj) , 1 + Mf+ (77))} , (4.1) 
then the first cigenpair of the Dirichlct problem (jl.lOp in the interval (0, 1) becomes 

Hi = 7r 2 , (ft (77) = 2~ 1/2 sin (ttt?) . 
Hence, the condition (|2. 15[) converts into 

> H± (i]) (cos 2 (7177) - sin 2 (ttt?)) dr\ = H± (77) cos {2irq) d-q. (4.2) 
Jo JO 

By Theorem [21 the two-dimensional problem (|2.3p ~ (|2.5p gets a trapped mode and the localization effect 
occurs in the case when the coefficient (|4.2[) in the Fourier series at least of one function H± € C [0, 1] 
is negative. 

For the symmetric thin domain (|4.1|) , i.e.. for H + = the simple first eigenvalue Ai = Af g 
ad (A±) of the operator A± of the problem (|2.3j) - (|2.5jl in II = H± must be treated as multiple in Theorem 
2] which gives the same asymptotic forms for A^ and A^ in p.7[) due to the equality A^ = A' 2 ' in (|3.3p . 
At the same time, the first eigenvalue A^ of the problem (| 1 . 3|) - (| 1 . 5[) is simple by the maximum principle. 
To display a difference in the asymptotics of A^ and A^, we need to construct the second term in the 
expansion of the eigenvalues. 

The Fourier method ensures that, for £ > cu, the eigenfunction U\ of the problem (|2.3p ~ (|2.5p admits 
the decomposition 

U x (£) = C x exp (- {tt 2 - Ax) 1/2 sin (ttt?) + O (exp (- (4tt 2 - Ai) 1/2 c) ) • (4.3) 

The factor C\ in (|4.3p does not vanish because in the case C\ = the next term 

C* 2 exp (- (4tt 2 - Ai) 1/2 C) sin (2^77) 

in the decomposition of U\ becomes main but changes sign. 

Following J37J §5-6], we accept the asymptotic ansatze for an eigenpair {\ h ,u h } of the problem 

A' 1 = h~ 2 (Ai + eA[ + ...) , (4.4) 



= ( 6 ±^i h ^ C 1 t ^))) + h- 1 (1 t ^) 



± 

where the dots stand for neglectible terms and 



e = exp [~2h- x (tt 2 - Aj) 1/2 c) . (4.5) 

We insert (I4.4p into the equation (| 1 . 3[) and the Neumann boundary conditions (jl.5p while remarking 
that the Dirichlct conditions (|1.4p are satisfied. In view of (|2.ip we write 

exp (- (tt 2 - A^ 1/2 C ± ) = exp (-h- 1 {it 2 Ai) 1/2 (1 T *)) = exp (-ft" 1 (tt 2 - A x ) 1/2 (2 - (1 T *))) = 



exp ( -2k- 1 {tt 2 - A 1 ) 1/2 ) exp ({tt 2 - A 1 ) 1/2 C T 
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Figure 8: Division of a symmetric domain 



that explains the form (|4. 3|) of the new small parameter. Furthermore, we derive the following problems 
in IT to determine the second terms in the asymptotic ansatze: 

-A f t/; (£) - AiC/; (£) = & ± A^ (0 , £ G n, (4.6) 
U' ± (0 = 0, € e S, 

a„C/i (0 - -M„ fcxp (V - Ax) 172 ^) sin (to?)) , f G r. 



Since the eigenvalues A' ± are simple, the only compatibility condition in the problem (|4.6| reads: 

A' ± 6± = A' ± 6± ^ C/i (0 2 dt = b T Ux (0 3 V (exp ((tt 2 - Ai) 1/2 C ± ) sin (ttt?)) rfr; =: 6 T F. (4.7) 

The authors do not know a way to confirm the assumption 

F^O. (4.8) 

However, under this assumption, the system of the two (±) linear algebraic equations (|4.7p has the 
following couple of solutions 

A' ± = —F, b± = ±2- 1 ' 2 and A' ± = F, b± = 2~^ 2 

that complete the asymptotic ansatze (|4.4[) . A straightforward and simple modification of the arguments 
in Section [3] justifies the constructed asymptotics and, in particular, gives the relations 

A? = h- 2 ( Ai - eF + O ( e 3 ' 2 ) ) and A* = h~ 2 ( A, + eF + O ( e 3 ' 2 ) ) . (4.9) 



The formulas (|4.9[) . (|4.8[) show the asymptotic fission of the first two entries in the eigenvalue sequence 
(|1.7p of the problem (|1.3p - (|1.5[) in the symmetric thin domain fl h . We note that Af and A§ , respectively, 
are eigenvalues of the spectral problem restricted on the half domain (see Fig. [5]) 

n h+ = {x = (y, z) : y G (0, h) , < z < 1 + H (h^y) } 

with the Neumann and Dirichlet conditions on the straight lateral side while, respectively, the even and 
odd in z extensions of the eigenfunctions complete eigenpairs of the problem (|1.6p in the domain (|4.1D . 
In this way one sees that the localization effect occurs at the both ends simultaneously. However, this 
observation does not help to verify the inequality 



4.2 The Dirichlet problem 

Let us consider the equation (|1.3[) with the Dirichlet boundary conditions on the whole boundary 

u h (x) = 0, x G dn h . (4.10) 
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Figure 9: Infinite cane- head domain 



By the classical approach [28], the Dirichlet problem 



- A e U = AU 



IL 



U = on dm 



(4.11) 



in the semi-infinite cylinder 



n± = {^iriew, C ± > -H±(ri)} 



(4.12) 



with the profile function H± G C 1 (uJ) , has no solution in H 1 (Ji± ; 911 ± ) . In other words, the unbounded 
operator A±, generated by the quadratic form (|2.7[) in H 1 {U±;dU±) (see [TH §10.2]), possesses empty 
point spectrum. As a result, we see that the localization effect discovered in the previous sections for 
the mixed boundary value problem (|1.3|) - (|1.5|l under the hypotheses of Theorems [2] and [3] does not occur 
in the same domain (|1.1|) . However, employing a result in [29], we observe the effect in the case of a 
dumbbell domain in Fig. [3[ 

Let G± C K™ = {£ : £ < 0} be a domain such that the boundary <9G± contains the set {£ : = 0, T] G o>} 
and the closure G± is compact. We set 



n± = G±U(wx [0,+oo)) , 

n h = (ux[-l,l])u\jG h ± , G£ = {s : £± e G±} 



(4.13) 
(4.14) 



(see Fig. [9]and[3l respectively). 

We assume that the first eigenvalue Af^_ of the Dirichlet problem in the bounded domain G± lies 
below the cut-off \i\ of the problem (|4.1ip . (|4.10p in the unbounded domain (|4. 12[) . One may fulfill this 
requirement by inflating a domain of a fixed shape. Following [29], we extend the first eigenfunction 
<&i by zero from G± onto Il-j- and apply the minimum principle ([TH Thm. 10.2.1]) to derive that 



min{A : A G o (A±) } 



inf 



a (U, U) 



UeHi{Tl ± ;dTl ± ) \\U]L 2 (II±) 



< 



^;i 2 (n±)| 



(n±)| 



(4.15) 



lkf ± ;i 2 (G±)| 



Af ± < ^1. 



As in Section [21 the relation (|4.15p assures that the discrete spectrum of operator A± is not empty 
and the problem (|4. 1 1|) in the cane-head domain (|4.13[) has an eigenvalue Ai G (0, ^i). In this way, by 
the inflation of G±, one can place any given number of eigenvalues in the interval (0, /ii). 

We repeat word by word the argumentation in Section [3] and arrive at the following assertion. 



Theorem 11 Let the operators A± of the problems j4-H\ ) in H± have the eigenvalues VS. 1}) . Then 
Theorem [7] keeps the validity for the eigenvalue sequence jj.7| ) of the Dirichlet problem il.3\) , \4-10\l in 
the dumbbell domain il h in ft4.14\ )- 
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Figure 10: Division of a symmetric domain 



4.3 The Neumann problem 

Let us change (|4.10| for the Neumann boundary condition 

d„u(x) = 0, x £ dfl h , (4.16) 

where d v stands for the differentiation along the outward normal v which is defined almost everywhere 
on the Lipschitz (by our assumption) boundary of the domain (|l.lj) (see Fig. [TJ. 
The integral identity 

(V x tt h , V x v) nh = X h (u'\v) nh , DGff 1 {n h ) , (4.17) 

posed an the whole Sobolev space H 1 (Q h ), serves for the Neumann problem (|1.3[) . (|4.10p . It is well- 
known that the eigenvalues 

= A£ < \'l < X% < ... < Xp < ... -» +00 (4.18) 
of the problem (|4. 17)) satisfy the asymptotic formula 



< ch^p 3 for h £ (0, hp' 4 ) , (4.19) 



where the positive constants c and h are independent of the eigenvalue number p £ N (see [25] §7.1] for 
the estimation of the bounds in (|4.19|) and also [30] [31] for similar results in other singular perturbed 
spectral problems). Note that M p = jir 2 p 2 , p = 0,1,..., are eigenvalues of the one-dimensional limit 
problem 

d 2 z w (z) = M p w {z),z£ (-1, 1) , d z w (±1) = 0, (4.20) 

obtained by the traditional procedure of the dimension reduction (see [3H Hi] and others) . 

For p = O (h~^ , formula (|4.19p no longer displays an asymptotics of the eigenvalue because the 

term jir 2 p 2 becomes of the same order h^ 1 as the bound chip 3 . As mentioned, e.g., in [Ml [3D], the 
above observation can help to reveal other asymptotically stable series of eigenvalues. 

To indicate such eigenvalues, we apply the method of artificial boundary conditions proposed in [35j 
to find out the point spectrum in the continuous spectrum of the Neumann Laplacian in a symmetric 
strip with a finite obstacle in the inside. 

Let us assume that the thin domain (jTTTJ) is symmetric with respect to the plane {x : y\ = 0} , i.e., 
formally we have 

Q h = {x : (-yi,V2,...,Vn-i,z) £ n h } . (4.21) 

Then we restrict the equation (jl.3[) on = {x £ Q h : yi > 0} and the boundary condition (|4.16[) on 
(dn h ) A = {x£ dn h : yi > 0} (see Fig. [IOJ while imposing the artificial boundary conditions 

u h A = on T' 1 = {x £ n h : yi = 0} . (4.22) 
The variational formulation of the new mixed boundary value problem refers to the integral identity 
(V^V^aW - A* «O n „ ,v £ H 1 {n h A ;T h ) . (4.23) 
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Let us make several observation. First, any cigenpair {A A ,7t A } of the problem (|4.23p becomes an 
eigenpair of the Neumann problem (|4. 17)) after the extension of u\ from f2^ onto VL h as an odd function 
in y\. Second, in the limit boundary value problem 

-A i U A = A A U A in n A = {£ G n : 771 > 0} , (4.24) 
C/ A =0on T A = {£ G n : m = 0} , d v U A = on 0II A \T* 

posed in the half of the cylinder (|2. 16|) . the continuous spectrum [/_t A , +00) begins with the first eigenvalue 
/i A > of the problem 

-A^ A = fi A ip A in lu a = {77 e uj : 771 > 0} , (4.25) 
ip A = on v A = {77 G uj : 171 = 0} , d„tp A = on &j a \JA 

(compare with ()2.3() - (|2.5[) and p. lip ). Third, owing to the positive cut-off /j a , the arguments used in 
Section |2~31 to prove Theorem [5] maintain the following assertion. 

Theorem 12 Let the cross-section uj be symmetric, i.e., ui = {?/ : (— 771, 772, ...,r) n -i) G uj}, and let the 
even in r/i function H £ C (uj) meet the condition 

H (v) (|V^ A (n)f - fi A ifi (r ? ) 2 ) d V < 0, (4.26) 

where (n A ,tp A ) is the first eigenpair of the problem \4-2h^ . Then the mixed boundary value problem 
\4-.24\j has an eigenvalue below the cut-off /i A for the continuous spectrum. 

Remark 13 We have used in {2.11^ the formula 

2 (V^x, <piVr,H) u = {<pt,\H) u (4.27) 

in order to derive the simplified condition V2.18)) in Theorem^ In general the change <px,uj 1— > (p A ,uj A 
makes fr4-27\ l wrong, since tp A does not vanish on v A . However, in the case V,,-ff (77) = for r/ G v the 
assertion in Theorem^ remains valid under the symmetry assumption £4 -21]) . I 

The Dirichlet conditions on the part v A of the boundary are sufficient to realize the justification 
scheme, developed in Section [3] Indeed, the key inequality (|3.9[) ought to be replaced by the inequality 

\\V v w;L 2 (uj a )\\ 2 > /j a \\w-L 2 (w A )f, w G H 1 (uj a ;v a ). 
Hence, the above observations provide the following assertion. 

Theorem 14 Let A± G (0,/i A ) be an eigenvalue of the problem \4- 6 M ) mlT± = {£ : 77 G ui A , £ > —H± (77)} 
(cf. Theorem \12\ and Remark \13]) . Then the Neumann problem \4-2Sfy has an eigenvalue A^^ such that 



\h 7,- 2 A A 



where Ca and r A are certain positive numbers. 



< c A exp (-/7 _1 r A ) (4.28) 



We emphasize that the eigenvalue number N (h) in Theorem [T?] depends on the small parameter 
h and N (h) — > +00, as h — > + . Moreover, we cannot assert that the only eigenvalue A^^ satisfies 
the relation (|4.28|) . These issues follow from the fact that the first series of eigenvalues with the stable 
asymptotics (|4.19[) is of order h a but the eigenvalue A^^ is of order h~ 2 . 
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